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Abstract We use observations of Sunyaev-Zel'dovich effect and X-ray surface brightness 
to reconstruct the radial profiles of gas temperature and density under the assumption of a 
spherically symmetric distribution of the gas. The method of reconstruction, first raised by 
Silk & White, depend directly on the observations of the Sunyaev-Zel'dovich effect and the 
X-ray surface brightness, without involving additional assumptions such as the equation of 
state of the gas or the conditions of hydrostatic equilibrium. We applied this method to the 
cluster RX J 1347.5-1 145, which has both the Sunyaev-Zel'dovich effect and X-ray observa- 
tions with relative high precision. It is shown that it will be an effective method to obtain the 
gas distribution in galaxy clusters. Statistical errors of the derived temperature and density 
profiles of gas were estimated according to the observational uncertainties. 
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1 INTRODUCTION 

Galaxy cluster is known to be the largest virial gravitational bound system in the universe. The strong 
gravitational potential heats the gas (Hydrogen and Helium) to be fully ionized with a very high temperature 
~ 10 7 — 10 8 K. The thermal electrons collide with ions and emit bremsstrahlung radiation, which is mainly in 
X-ray band, making galaxy clusters strong X-ray sources. The intensity of thermal bremsstrahlung radiation 
is approximately proportional to Te^ 2 n 2 JBirkinshaw 19991 1 . So by measuring the X-ray surface brightness, 
we can obtain information about the gas temperature and density disnibutions of the galaxy cluster. One 
of the most popular models is the so-called isothermal (3 model (Cavalier e & Fusco-Femiano 1976b , in 
which the temperature disnibution as a function of radius is assumed to be constant and the density profile 
takes the form n e (R) = n e o(l + -R 2 /i? 2 )~ 3/3 / 2 . This simple model gives good approximations to many 
galaxy clusters. However, more and more observations indicate that the temperature distribution may not 
be a constant over the whole galaxy. Hughes et al. (1988) introduced a revised model with a constant 
temperature in a range r < r; so and a decrease at large radius, while the electron density has a truncation 
at a radius riim. This form provides an adequate description for some of the clusters. Further modification 
on the (3 model has been proposed to deal with the existence of cooling flow in the center of the cluster 
(IFabian et al. 19841 >. 

In order to get rid of the degeneracy between T e and n e , further assumptions like the hydrostatic 
equilibrium condition or the poly tropic equation of state (EOS) of the electron gas dCowie et al. 19871 
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IXue & W u 2000, Wu & Chiueh 200Tb . or the spectrometric analysis ( IJia et al. 20061 1 are needed. However, 
the results are still model dependent. A realistic reconstruction is expected to come directly from the obser- 
vations. 

The inverse Compton scattering between the electrons in clusters and the cosmic microwave 
background (CMB) photons, namely the Sunyaev-Zel'dovich (SZ) effect ( |Zeldovich & Sunyaev 1969} 



Sunyaev & Zeldovich 1970a b), provides a possible way to achieve this goal. The distortion of the CMB 
spectrum is related to the electron temperature and density, and will provide another equation of T e and n e . 
Combining the SZ effect and the X-ray observations, one can obtain the electron temperature and density 
profiles. This method was first suggested by Silk & White (1978)| but has not yet been put into practice due 
to the limited instrumental sensitivity and resolution in the detection of the temperature fluctuation of CMB. 
Up to the 1990s some preliminary observational results of the radial temperature distribution were obtained 
for a few clusters (Birkinshaw et al. 199U Birkinshaw & Hughes 1994 1, but it was still difficult to obtain 



precise reconstructions of the temperature and density profiles because the data points were too sparse and 
uncertain dWu 20031 ). 

In this paper, we try with this method to acquire some preliminary information about the gas temper- 
ature and density distribution. For simplicity, we assume a spherically symmetric distribution of the gas in 
the clusters. We adopt the ACDM model with Q M = 0.27, = 0.73, and Hubble constant H = 71 km 
s _1 Mpc -1 (Spergel et al. 2003 1. This paper is organized as follows: we describe the method of reconstruc- 
tion in Section|2l In Section |3J we apply this method to the galaxy cluster RX J 1347.5-1 145 to reconstruct 
its temperature and density profiles, together with an estimate of the uncertainties. The conslusions and 
discussions are presented in Section|4] 

2 METHOD OF RECONSTRUCTION 

The fluctuation of the temperature of CMB due to the thermal SZ effect is (Birkins haw 19991 ) 

— ^ =9{x)y{R), (l) 

where g(x) = x coth(x/2) — 4 with x = hv IksTcMB the dimensionless frequency, y(R) is the 
Comptonization parameter as a function of the two-dimensional skymap radius R. In the Rayleigh- Jeans 
limit, i< 1, g(x) = —2, so we have ATcmb(R)/Tcmb = — 2j/(i?). The Comptonization parameter is 
related to the electron temperature and density by 

f°° rdr 
y(R) = 2A y / T e (r)n e (r)^==, (2) 
Jr vr z -R z 

where A y = kB(JT/m e c 2 , fcs is the Boltzmann constant, ot the Thomson cross section, m e the mass 
of electron and c the speed of light, r denotes the three-dimensional spatial radius from the center of the 
galaxy cluster, while R = cIaO with cLa the angular diameter distance from the Earth to the cluster and the 
(projected) angular separation. T e (r) and n e (r) represent the electron gas temperature and number density 
as functions of the radius r, respectively. 

The X-ray surface brightness of thermal bremsstrahlung emission is 

S *W = a n I H ' 2A * / Tl' 2 {r)nl[r)C{T e ) (3) 
47r(l + z) 4 J R y/r 2 - R 2 

where A x = 3 j£f c s ( 3 2 ^ fc g; ) 1/2 ^ e g, Me = 2/(1 + X) with X = 0.768 the primordial Hydrogen mass 
fraction, g w 1.2 the average Gaunt factor (this average leads to an uncertainty of ~ 20% in the bolometric 
emissivity, see e.g., IEttori 20001 ) and z the redshift of the cluster. C(T e ) = exp[— -B m i n (l + z)/fcsT e ] — 
exp[— £ ; max (l + z)/ksT e ] is a correction factor from the whole band of thermal bremsstrahlung emission 
to the detector bands [E min , E max ] (e.g., for ROSAT, E min = 0.1 keV and S max = 2.4 keV). 

Using Abel's integral equation, Equations © and @ can be inverted to (Yoshikaw a & Suto 19991 ) 

r . (rK(r) _ » (4 > 

IT Ay Jryz AR y/ R 2 - T 

^ r r dS x (R) dR 

7T- (P) 



T^ 2 (r)n 2 e (r)C(T e ) = *<L 



a x dR ^/m 
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Table 1 Parameters used 
and O. 


in Equations 




X-ray S^oCergs s _1 cm -2 arcmin -2 ) 
1.34 x 10" 10q 


# ca: (arcsec) 
4.29 ±0.10 6 


0.535 ± 0.003 


SZ yo 

4.1±£;J x 10- 4 


9 C y(arcsec) 
56.31™ 


(3cy 

0.89«; 2 ? 



"Derived from the Chandra AC IS-S detector counts, see 
http://heasarc.nasa.gov/Tools/w3pimms.html 

"Note that the angular diameter distance is different from Allen et al. (2002)] 
due to different cosmological models. 



Thus, given the observational distributions of y(R) and S X (R), we can easily find the electron temperature 
and number density distributions from Equations and (0. 

3 APPLICATION TO CLUSTER RX J1347.5-1145 

RX J1347.5— 1 145 is a highly X-ray luminous and dynamically relaxed galaxy cluster with redshift 
z = 0.451. Both ROSAT and Chandra have measured the X-ray emission of this cluster with high pre- 
cision (Schin dler et al. 19971 Allen et al. 2002). The X-ray surface brightness profile can be parameterized 
following the conventional (3 model as 

*W = *(l + ^) , (6) 

where dA is the angular diameter distance of the cluster. For our adopted cosmological model, we find 
dA = 1185 Mpc. Using the 18.9 ks Chandra ACIS exposure, |Allen et al. (2002)| gave a detailed X-ray 
image of this cluster in the energy band 0.3 — 7.0 keV. After subtracting the south-east excess, the X-ray 
surface brightness can be well described by Equation ((6). 

Detections of the SZ effect at different frequencies in this cluster were published by several groups 
(Pointecouteau et al. 2001; [Komatsu et al. 20011IReese et al. 20021 >. An empirical formula similar to the (3 
model can also be used to parameterize the observational Comptonization parameter of the SZ effect y(R), 

We use the fitting parameters from |Pointecouteau et al. (200 1)| The parameters in Equations (O and (0 are 
listed in Table[3] 

Substituting y(R) and S X (R) in Equations (01 and © by Equations © and (0, we can obtain 



T e (r)n e (r) 



1 y m v /2 - 1/2) T(3(3 cy /2) 



-3/3c„/2 



^A v d A e cy T(Wc V /2 + 1/2) V d\6 2 cy 
= C y , (8) 

T^(r)nj(r)C(T e ) = V^ 1 + ^ S MW CX 1/2) T(3p cx ) 



A x d A e cx r(3/3 ca + 1/2) V d? A 0l x , 

= C x . (9) 

Then 

T e {rfl 2 /C(T e ) = C 2 y /C x , (10) 
n e (r) = C y /T e (r). (11) 

From Equation ( fT0T > we know that, if 6 CX = 9 cy , /3 CX = (3 cy , the temperature should be independent of 
the radius, as expected from the isothermal f3 model. If (3 cy > (3 CX , the temperature will decrease at large 
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radius; while for (3 cy < f3 cx , the temperature will increase in the outer region of the cluster, which seems to 
be unreasonable. 

It is easy to calculate T e (r) and n e (r) using the parameters given in Table|3] To obtain the uncertainties 
of the derived temperature and density profiles, we run a Monte-Carlo (MC) sampling of the parameters 
according to their uncertainties. The parameter is thought to be Gaussian distributed peaked at the central 
value with width of ler error-bar. However, we note that the ler error-bars are not the same in the "+" and 
"— " directions. So the distribution of a parameter is the combination of two Gaussian functions, connected 
at the peak point. For example, for a parameter a = ao-^' tne width of a > ao is o\ and a < ao is 1x2, 
with the connecting condition fci/eri = ki/a2- The probability distribution function of parameter a can be 
written as 

fci / (a— an) 2 \ . 

, ^ ex P(~ S^f") a>a °< M ~ 
^exp(-^^) a<a . 

Here, fci and &2 can be derived according to the normalization condition J p(a)da = 1. Some cut conditions 
are adopted in the MC sample. Equations (O and (0 require (3 cy > 1/3 and f3 cx > 1/6; the other parameters 
are required to be greater than zero. Even so, some of the parameter combinations make Equation dTOb 
unresolved. These cases were also discarded. 




10 1 10 2 10 3 

r/kpc 



Fig. 1 Derived profiles of temperature (upper panel) and density 

(lower panel) together with their uncertainties. Also shown are the results of some previous studies, see the text for 

explanations. 
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The average values and variances of the temperature and density of the MC sampling are shown in 
Figure Q] Shaded regions are the lx and 2x sample variances as represented by the ±ler and ±2<r er- 
rors. The results of previous authors are also shown for comparison. The solid crosses are the deprojected 
temperature profiles from the Chandra X-ray spectroscopy (Allen et al. 2002), the dotted crosses are 
the derived results from the SZ effect and X-ray data using a method slightly different from this work 
( |Kitayama et al. 2004 1, and the dot-dashed line is the extrapolated result from the Chandra spectroscopy 



observations (Schmidt et al. 2004). It is shown in this figure that our derived results are roughly consistent 
with the other results at 2a level. For the inner region of the cluster, our results show an under-estimation of 
the temperature. This might be due to the systematic errors of the observations (especially the SZ measure- 
ments). From Equations ( [Tol l and (flTT i we know that a larger yo or a smaller S x o may result in a larger central 
temperature. From |Pointecouteau et al. (2001)| it was shown that different fitting models could indeed give 
different yo- The density profile, as shown in the lower panel of Figure Q] is consistent with the results in 
|Schmidt et al. (2004)| for r < 500 kpc. As pointed out in |Allen et al. (200"2)] the surface brightness showed 
a steepening with increasing radius, which means an increase of the f3 parameter. A broken power law with 
(3 = 0.54 changing to 0.78 at r = 487 kpc could well describe the observations. We adopt a uniform (3 
parameter here, so for large radii our results are somewhat high. According to Equation (fTOl . a larger j3 cx 
will also lead to a higher temperature at large radius. 

4 CONCLUSIONS AND DISCUSSION 

Using the observational data of the SZ effect and X-ray surface brightness, we applied the method suggested 
by | Silk & Whit e (1978) to the galaxy cluster RX J1347.5-1 145, to reconstruct its gas temperature and den- 
sity profiles. This is a direct way to obtain the temperature and density profiles of galaxy clusters, without 
additional theoretical assumptions. However, the quality of the observational data strongly affect the recon- 
struction results. Our attempt on cluster RX J1347.5-1 145, which has both the X-ray and SZ measurements 
with relative high precision, demonstrates the effectiveness of this method. The derived results show similar 
behaviors as the previous studies. It indicates that there is a cooling flow in the center of the cluster, but the 
central value we derived is lower than the others. Poor quality and possible systematic uncertainties of the 
SZ effect data might be responsible for this discrepancy. It should be noted that the model parameters of the 
X-ray surface brightness and SZ effect were derived from finite area images around the center of the cluster, 
then extrapolated to large radii. This may result in additional uncertainty (see the discussion in Section[3]l. 
When high quality measurements of the SZ effect become available in the future, it will be a powerful tool 
to study the intrinsic gas distribution independent of any theoretical arguments. 
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